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ABSTRACT 

Despite a surface dominated by carbon-based life, the bulk composition of the 
Earth is dramatically carbon poor when compared to the material available at 
formation. Bulk carbon deficiency extends into the asteroid belt representing a 
fossil record of the conditions under which planets are born. The initial steps of 
planet formation involve the growth of primitive sub-micron silicate and carbon 
grains in the Solar Nebula. We present a solution wherein primordial carbon 
grains are preferentially destroyed by oxygen atoms ignited by heating due to 
stellar accretion at radii < 5 AU. This solution can account for the bulk carbon 
deficiency in the Earth and meteorites, the compositional gradient within the 
asteroid belt, and for growing evidence for similar carbon deficiency in rocks 
surrounding other stars. 

Subject headings: ISM: abundances - Stars: formation - Physical Data and 
Processes: astrochemistry 

1. Introduction 

The early stages of planet formation rely on the growth of small grains into larger bodies. 
Locally, the relative chemical composition of planets, comets and asteroids within the solar 
system provides indirect clues to the environment and processes that existed and operated 
billions of years ago. In this regard, the bulk elemental composition represents a key fossil 
record, as it can be compared to the distribution of elements available at birth as represented 
by the Sun, interstellar clouds, or extra-solar protoplanetary disks. 
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Figure 1 illustrates the dramatic gradient in the amount of carbon relative to silicon 
incorporated into planets, comets, and asteroids, represented by meteorites. In this figure 
the H/Si ratio serves to separate the various bodies with the proviso that the Sun/ISM 
represents the total amount of material available at the start. Material that formed far from 
the Sun (comets) incorporate water and other hydrogen rich ices, but do not include the H2 
gas that is the major component of the Sun and ISM. In contrast rocks within the terrestrial 
planet region (r < 5 AU) incorporate successively less H-rich ices, a well known result due 
to the snow- line lying inside the asteroid belt. Icy bodies have an abundance of carbon in 
the form of ices and amorphous carbon grains. In contrast the Earth has three orders of 
magnitude less carbon (relative to Si) than was available during formation as represented 
by the Sun/ISM. While the Earth measurement docs have uncertainty, as we do not know 
the amount of carbon sequestered in the core, estimates of the amount of hidden carbon still 
lead to significant carbon deficit (Allegre et al. 2001). Moreover, primitive carbonaceous 
meteorites, as tracers of the starting materials for the Earth, also exhibit carbon deficiency. 

In ISM ~ 60% of the cosmic carbon is locked in some form of graphite or amorphous 
carbon grains (Savage & Sembach 1996). These primordial grains need to be destroyed 
prior to planet formation to account for the composition of the Earth and carbonaceous 
chondrites. That carbon grains have been detected in comets with near-solar abundances 
(Fig. 1) suggests that any destruction mechanism was inactive in the cometary formation 
zone near the giant planets. Another hint lies in the variation of spectral classes within the 
asteroid belt, which may be indicative of compositional gradients in the relative amounts of 
carbonaceous material (Gradie & Tedesco 1982). Thus there are hints in the solar system of 
processes that were active in the hot inner nebula that kept primordial carbon grains from 
being incorporated into rocks. However, the same process was inactive in the colder outer 
nebula. 

In addition, the detection of carbon and iron in the atmosphere of some white dwarfs 
(Jura 2006) is an unexpected result as the strong gravitational field of a white dwarf is 
anticipated to lead to settling of heavy elements in a photosphere dominated by either 
hydrogen or helium (Paquette et al. 1986). The measured C/Fc abundance ratio in the 
atmosphere of these white dwarfs is well below solar and consistent with external supply by 
rocky bodies that are deficient in carbon (Jura 2008). In sum, the solution to the carbon 
deficit in the inner solar system is a natural outcome of star and planet formation, perhaps 
throughout the Galaxy. 

The standard model for the Earth's composition is that the planet formed from material 
less volatile than sodium, but with an additional contribution of volatiles (Waenke & Dreibus 
1988). In this approach, carbon is considered a volatile element since, in thermodynamic 
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equilibrium, in a mixture with a solar composition, it forms molecules such as CO and CH4. 
Therefore, the standard cosmochemical model presumes that the inner solar system was 
entirely vaporized at a temperature greater than 1,800 K (Grossman 1972). However, the 
computed disk midplane temperature at 1 AU is less than 1,000 K in most models of accreting 
pre- main- sequence stars (D'Alessio et al. 2005). Thus, primordial carbon grains should have 
been incorporated into the Earth. Carbon must have been present in the nebular gas in a 
volatile form (CO, CH4, or CO2). At 1 AU this gas would not be assimilated into silicate-rich 
planetesimals, but would be accreted by the Sun and gas giants or lost to space. Beyond the 
nebular snow-line large bodies could incorporate the carbon grains and ices, which can be 
contributed to the young Earth by impacts. Some previous solutions to this problem include 
ultraviolet heating and vaporization of infalling grains (Simonelli et al. 1997), reactions with 
OH (Gail 2002), and impacts (Bond et al. 2009). However, the infall phase is short and 
models suggest that either none or all grains, including silicates, are vaporized (Simonelli et 
al. 1997) and in the warm inner disk OH may not be abundant enough (Glassgold et al. 
2004; Glassgold et al. 2009; Woitke et al. 2009; Gorti & HoUenbach 2004). 

In this letter we meld measurements of relative amounts of carbon and silicon found 
within solar system bodies with our current understanding of the physics and chemistry in 
extra-solar protoplanetary systems to demonstrate that the innermost regions of the Solar 
Nebula was on fire and rocks forming within the terrestrial planet zone for the Sun and, 
likely other stars, are carbon-poor. 



The paradigm of terrestrial planet formation involves the growth of solids in the solar 
nebula, which is a protoplanetary disk. Grains initially present in the disk upper layers settle 
to a dust-rich midplane where they coagulate and grow into planetesimals (WeidenschiUing 
& Cuzzi 1993). Depending upon their size, grains can be suspended in the disk atmosphere 
via interactions with the turbulent gas (WeidenschiUing & Cuzzi 1993). In this theory sub- 
micron sized grains are held in the upper layers, by turbulent stirring, and coagulate to larger 
sizes where they eventually settle to the midplane. Grains settle to a vertical layer, Zgett, 
that depends on the disk physical structure, the strength of turbulence, parameterized by 
an a viscosity (Hartmann 2001), and the grain mass to surface area ratio. This height can 
be determined by setting the timescale of settling, tgett, equal to the timescale of diffusive 
turbulent stirring of the disk, tsUr (DuUemond & Dominik 2004): 
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Model Calculations and Results 
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Sc is the Schmidt number which relates the decouphng of grains from turbulence, m/a is 
the grain mass-to-cross section ratio, which is m/a = Ap^a/S, where pd is the grain density 
and a is the grain radius. Here, pd is calculated from graphite, whose density is 2.24 g cm~^, 
with assumed porosities (Ormel et al. 2007). Hp is the pressure scale height which is given 
by Hp = Csfl/^k- = \/G~Mj^7^ denotes the local Kepler frequency and the sound 
speed at the midplane, Cgfi = \J'^^-- For perfectly coupled particles S'c = 1, which is the 
case for grains smaller than a = 100 pm (DuUemond & Dominik 2004). When ^ = 1 the 
settling layer {zgett), where grains are suspended by turbulence, is reached. The depletion 
layer (zdepi), above which the disk is entirely depleted (by settling) of grains of that size, 
can be determined by setting ^ = 100 (Dullemond & Dominik 2004). We adopt a = 0.01, 
which seems appropriate according to the Balbus-Hawley instability mechanism to generate 
viscosity (Hartmann 2001). 

The gas and dust temperature and density in the disk are calculated by using the 
parameters listed in the caption of Figure 2 (see Nomura et al. 2007 for details) Q For a 
dust model, we use a dust size distributions of dn/da oc a~^^^ with grain radii ranging from 
0.01 pm to 10 pm, which are larger than those of the ISM grain model (Mathis et al. 1977), 
assuming the possibility of grain growth in dark clouds. In the calculation of stirring depths, 
we adopt various porosities for different grain sizes according to Ormel et al. (2007). 

Figure 2 presents the estimates of Zgett and Zdepi and the gas and dust temperature 
profiles at 1 AU in our disk model. We also show the chemical abundance profiles of key 
oxygen reservoir^. Disk accretion onto the star and an active chromosphere power X-ray 
and ultraviolet radiation fields that heat the disk surface to high temperatures at the depth 
where grains are suspended by turbulent stirring (Glassgold et al. 2004; Nomura et al. 2007); 
in contrast the disk midplane is cooler. A key facet of disk thermal models is that gas and 
dust temperatures are decoupled in the upper layers. Figure 2 shows that Tgas between Zs^tt 
and Zdepi can be quite high (> 500K) for grains < 10pm. In addition, within this depth range 
oxygen-bearing molecules are predominantly photodissociated with no = 1.8 x 10~^n// at < 



^The gas heating via PAHs is included in the model (e.g., Aikawa & Nomura 2006), but we note that the 
temperature structure may have some dependency on models (e.g., Gorti & HoUenbach 2004). 

^The chemical/thermal model did not include the chemistry of water and OH formation. For these 
abundances we used the analytical approximations of Bethell & Bergin (2009), which rely on knowledge of 
the local radiation field and gas temperature and density. The result is consistent with other self-consistent 
model results (e.g. Gorti & Hollenbach 2004). 
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5 AU. At 10 and 20 AU, however, the dust temperature is lower than the water desorption 
temperature (110 K) at all depths so that a significant amount of oxygen is locked on grain 
surfaces, as expected beyond the snowline. 

Hot oxygen atoms erode carbon grains and release the carbon into the gas. The physical 
and chemical conditions between Zsett and Zdepi at the terrestrial planet forming zone in our 
model are such that the chemical sputtering of carbon grains by hot oxygen atoms will be 
operative. Thus, for carbon, the need to increase the dust temperature above the sublimation 
point, as in the standard cosmochemical model, is evaded. Silicates will not be subject to 
similar chemisputtering (Barlow & Silk 1977) and remain the seeds for terrestrial worlds. 

The oxidation rate of carbon-based grains is: koxy — n-oo^voY (s~^). uq and Vq arc the 
number density and the thermal velocity of oxygen atoms, respectively. Y is the yield (carbon 
atoms removed per incident oxygen atoms). We adopt laboratory studies for the yield as 
given in Table 2 of Draine (1979), which covers the appropriate temperature range. For 
Tgas — 500 K and 1000 K, y ~ 0.02 and 0.2, respectively. A summary of recent experiments 
at Tgas > 1000 K can be found in Vierbaum and Roth (2002); these results are consistent 
with our assumptions. 

The key question is the destruction timescale, tdest, as grains cascade to smaller sizes 
when carbon atoms are removed. If tdest is shorter than tsur, any carbon grain that is lifted 
above the settling layer will be destroyed. To determine tdest we have constructed a kinetic 
chemical network, such as + 0^gn + C + (with n = # of carbon atoms) with the 
rate specified above. This network is then extended down to carbon atoms and we solve 
the system of ordinary differential equations to calculate tdest of the biggest grain (0.002 
yum) feasible to the computational memory and time. We then extend these results to larger 
grains by extrapolating this timescale (see Eq. (4) below). 

As shown above, the reaction rate (koxy) is proportional to the grain surface area {a = 
TTO^) in a given temperature and abundance of atomic oxygen. For a, we have related the 
grain volume to its mass: 



i-' = ^. (3) 

<J Pd 

Here mc is the total mass of carbon atoms contained in the given size grain. Based on this 
model with no porosity, a 0.002 fim grain {g?,ooo) contains 3000 carbon atoms. For the initial 
abundance of each grain (g^) in this kinetic calculation, we adopted the total number of 
carbon atoms (per H atom) derived from the MRN interstellar grain size distribution for 
graphite (see Figure 6.19 of Tielens 2005) and extrapolated to very small grains. 
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Figure 3 shows the abundance variation of atoms in carbon grains as a function of time 
at 1 AU. The physical conditions for the kinetic calculation are obtained from our disk model 
at the depth (0.122 AU) where Tgas = 500 K. These results demonstrate that, assuming no 
source term, there is a rapid decay of carbon grains to reduced carbon content and size. 
We define tdest of a certain size of grain when X{gi) decreases by 4 orders of magnitude in 
the reservoir. Therefore, 0.002 /im grain ((73000) is destroyed within 2 years at the depth of 
Tgas = 500 K at 1 AU. tdest of 0.002 /xm grain at 3, 10, and 20 AU are 14, 100, and 400 
years, respectively (see Table 1). In the calculation of tdest at 10 and 20 AU, we assume all 
oxygen is atomic although a significant amount is likely locked on grain surfaces as water 
ice. Therefore, the calculated t^est at 10 and 20 AU is the lower limit. In this disk model, 
the same size grains are destroyed at the depth where Tgas = 1000 K about 3 times faster 
compared to the depth where Tgas = 500 K. Our estimates of destruction timescales are very 
conservative as we have examined the timescale at 500 K and defined the timescale after a 
several orders of magnitude decrease in abundance. 

As grains grow in size by colliding aggregation, the porosity, the measurement of volume 
that is not occupied by carbon atoms, increases (Ormel et al. 2007). The enlargement 
parameter ip is defined hj ip = V/V*, where V* is the volume that carbon grains occupy 
in its compacted state while V is the total spatial extent of the given size grain. Thus, the 
porosity is defined as[{V — V*)/V]xlOO — [1 — xlOO. Therefore, the same size of grain 
of a higher porosity contains a lower number of carbon atoms to reduce tdest- ip —1 does 
not necessarily mean no porosity, but we assume no porosity for =1 for the conservative 
calculation of tdest (i-©-) the longest destruction timescale) of a given size grain. According 
to the results of Ormel et al. (2007), we adopt ip =1, 8, and 27 for 0.1, 1, and 10 /xm 
grains, respectively. Therefore, the porosity of 0.1, 1, and 10 /im grains are 0, 88, and 96%, 
respectively. 

tdest at a given condition is proportional to 1/koxy x H^{C)graim where ^{C)grain is the 
total number of carbon atoms contained in the grain, and thus, proportional to the size of 
grain (l/a x oc a). However, the porosity reduces tdest since 4i^{C)grain of a given size 
decreases. Therefore, tdest of a given size grain can be calculated by: 



tdest,a — (q um^^" ^<iest,0.002, (4) 

where tdestfi.om and tdest,a are tdest of 0.002 /xm grain and a grain with size a. ipa is the 
enlargement factor for grains of a certain size. The results of the combination of grain 
size and porosity are presented in Table 2. At the depth of Tgas = 500 K, tdest < tstir 
for grains smaller than 10 //m at 1 and 3 AU while at 10 and 20 AU, only very small 
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grains (0.01 /im) satisfy the condition. Therefore, all primordial carbon grains are destroyed 
within the terrestrial planet forming zone by turbulent stirring above Zgett within tsur- This 
effectively stops incorporating carbon grains in the initial stages of planetcsimal formation 
- the coagulation of small grains in the upper layers and the subsequent settling to lower 
heights. Moreover, this mechanism will decrease in efficiency within the asteroid belt, which 
could produce the observed asteroid compositional gradients. 



3. Discussion 

Carbon grain destruction in the dense UV free midplane has been explored by Finocchi 
et al. 1997, Bauer ct al. 1997, and Gail (2001, 2002). In this work OH is the primary agent 
for carbon grain destruction. However, the amount of OH in the disk midplane has some 
uncertainty as these authors determine an OH abundance of ~ 10"^'' — 10~^ relative to H2 
(Finocchi et al. 1997, Bauer et al. 1997), while Glassgold et al. (2009) find an abundance 
well below 10~^^. Our model refies on an abundance of oxygen atoms on the disk surface 
(e.g. Fig. 2) in a region that is clearly exposed to UV radiation, ensuring a ready supply of 
atoms with 0/OH ^ 1. However, our model does not preclude destruction of carbon grains 
in the midplane by OH, if the abundance is high enough. 

As a result of chemical destruction of carbon grains all carbon in the inner solar system 
will reside in the gas, which will accrete into the star or dissipate. Carbon grains that reach 
sizes of a > 50 — 100 /im in the outer disk will settle below the layer where chemisputtering 
dominates. This brings up an important question: if the Earth did not incorporate carbon 
grains or accrete sufficient amounts from the gas, how did the Earth assimilate any carbon? 
We suggest that maybe the problem of the Earth's carbon is related to the origin of the 
Earth's oceans. At 1 AU it is likely that grains were too hot to be sufficiently coated with 
ices and instead the Earth potentially received its water from impacts of bodies that originate 
beyond the nebular snow-line near 2.5 AU (e.g. Morbidelli et al. 2000, but see also Drake 
2004). These bodies are hkely also carbon rich (e.g. Gradie & Tedesco 1982) and could have 
brought carbon with water ice. The mass of carbon in the Earth's mantle is ~ 5.5 x 10^^ g 
with perhaps up to an order of magnitude additional carbon sequestered in the core (Allegre 
et al. 2001). In the study of terrestrial planet formation performed by O'Brien et al. (2006) 
they estimate that ~ 15% (median value of simulations) of the Earth was delivered by bodies 
that originated beyond 2.5 AU. If we assume that these bodies contained the same amount 
of carbon as in carbonaceous chondrites of ~ 10 mg/g (Waenke & Dreibus 1988) then this 
corresponds to a delivery of roughly 10^^ g of carbon, which is in the exact range of that 
estimated for the Earth (see Bond et al. 2009 for an additional look at this issue). 
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Fig. 1. — Carbon abundances relative to Silicon in the Solar System. Abundance ratios 
of carbon (listed as C) relative to silicon for representative bodies in the solar system as 
a function of their H/Si content. Shaded areas denote specific types of bodies (2 Oort 
cloud comets, 4 classes of chondritic meteorites, Earth) or instances with strong similarities 
(Sun and ISM). Abundance estimates are taken from the literature (Geiss 1987, Wasson & 
Kallemeyn 1988, Sofia et al. 1994, Allegre et al. 2001, Min et al. 2005, Grevesse et al. 
2007). The amount of hydrogen on the Earth is estimated assuming a hydrosphere mass of 
1.4 X 10^*^ g. Figure adapted and expanded from Geiss 1987. 
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Fig. 2. — Left: Temperatures and Zgett and z^epi in the protoplanetary disk atmosphere. 
Gas and dust temperatures as a function of vertical height above the midplane at 1 AU. 
Also shown are the estimated depletion and settling depths for various grain sizes. Right: 
Chemical abundances of key oxygen carriers as a function of vertical depth. The chemical 
calculation is stopped at z ~ 0.1 AU when the gas/dust temperatures reach equilibrium. 
The model assumes an axisymmetric a-viscous disk (e.g., Hartmann 2001) with a = 0.01 
and M = lO^^M© yr~^, which is surrounding a central star with parameters of = I.OM0, 
i?* = 2Rq, and T^, = 4000K. Both UV and X-ray spectral profiles from the central star are 
assumed to be the same as those of TW Hydra, and they are normalized so that the total 
FUV luminosity is Lpuv = 1 x 10^^ ergs s^^ (which corresponds to Gpuv = 300 Go at the 
disk radius of 100 AU, where Go = 1.6 x 10~^ ergs s~^ cm~^) and the X-ray luminosity is 
Lx = 1 X 10^^ ergs s~^ 
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Fig. 3. — Kinetics of carbon grain destruction at 1 AU. Plot of grain abundances - maximum 
grain size (a(g'3ooo)) is ~ 0.002 /xm - with variable number of carbon atoms as a function 
of time using conditions at the depth where Tgas = 500 K (see Table 1). Various line types 
represent grains containing different carbon atoms. The gray arrow indicates the timescale 
{tdest of 0.002 iim grain) in which the abundance of gi decreases by 4 orders of magnitude. 
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Table 1: Destruction timescale of 20 A grains at the depth of Tgas — 500 K 





1 AU 


3 AU 


10 AU 


20 AU 


Depth (AU) 


0.122 


0.494 


2.047 


5.180 


Hp (AU) 


0.033 


0.101 


0.460 


1.164 


n(H2) (cm~^) 


1.42 X 10^ 


2.03 X 10® 


2.50 X 10^ 


7.50 X 10^ 


tstir (yr) 


160 


840 


3200 


8400 


tdest (yr) 


2 


14 


100 


400 



Table 2: Destruction timescale of various size of grains at the depth of Tgas — 500 K 



Grain Size (//m) 


Porosity (%) 


1 AU 


3 AU 


10 AU 


20 AU 






(160 yr)^ 


(840 yr)^ 


(3200 yr)^ 


(8400 yr)^ 


0.01 





10 yr 


70 yr 


500 yr 


2000 yr 


0.1 





100 yr 


700 yr 


5000 yr 


20000 yr 


1 


88 


125 yr 


875 yr 


6250 yr 


25000 yr 


10 


96 


370 yr 


2593 yr 


18520 yr 


74000 yr 



"'tstir at each radius as given in Table 1. 



